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ABSTRACT 


It  is  well  known  that  physical  processes  occurring  at  the  air-sea  interface,  such  as  air-sea  mass, 
momentum  and  energy  exchanges,  are  controlled  by  the  ocean  surface  roughness  that  has  often  been 
characterized  in  the  literature  as  being  homogeneous  in  space  and  time  from  a  large-scale  point  of 
view.  However,  the  ocean  surface  roughness  is  heterogeneous  within  a  wavelength.  The  surface  areas 
on  wave  crests  are  usually  rougher  than  those  on  wave  troughs.  The  mechanisms  range  from  capi  ary 
ripples  on  wave  crests  to  larger  scale  wave  breaking.  Ambient  swell  and  non-uniform  currents  can 
also  cause  significant  spatial  variations  of  the  surface  roughness  as  short  waves  become  steeper  on  the 
swell  crest  than  those  on  the  swell  trough.  The  study  parameterizes  the  spatial  variability  of  the  sea 
surface  roughness  and  quantifies  the  resultant  variations  of  momentum  and  energy  exchanges  between 
the  atmosphere  and  the  coastal  ocean  using  a  phase-resolving  numerical  wave  model.  Numerical 
experiments  have  been  carried  out  to  test  a  hypothesis  for  the  suppression  of  the  wave  growth  rate  by 
the  following  swell.  A  simple  mechanism  is  suggested  to  explain  the  influence  of  swell  on  the  short 
waves  observed  in  the  numerical  and  laboratory  experiments. 
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1.0  INTRODUCTION 

Many  naval  operations  occur  over  and  in  the  atmosphere-ocean  interface.  The  success  of  such 
operations  depends  on  the  understanding  of  the  structure  and  behavior  of  the  air-sea  interface.  It  is 
well  known  that  many  physical  processes  occurring  at  the  air-sea  interface,  such  as  aii-sea  mass, 
momentum  and  energy  exchanges,  are  controlled  by  the  ocean  surface  roughness  that  has  often  been 
characterized  in  the  literature  as  being  homogeneous  in  space  and  time  from  a  large-scale  point  of 
view.  Recently,  it  is  realized,  however,  that  the  utilization  of  remote  sensing  techniques  for  the  study 
of  the  air-sea  interface  requires  the  in-depth  understanding  of  the  spatial  and  temporal  variability  of 
the  ocean  surface  roughness  at  multiple  scales. 

The  ocean  surface  roughness  is  indeed  inhomogeneous  within  a  wavelength.  The  surface  areas  on 
wave  crests  are  usually  rougher  than  those  on  wave  troughs.  The  mechanisms  range  from  capillary 
ripples  on  wave  crests  to  larger  scale  wave  breaking  (see  e.g.,  Longuet-Higgins,  1963;  Jiang  et  al., 
1 999;  Ducan  et  al.,  1 999).  An  ambient  swell  or  current  can  also  cause  significant  spatial  variations  of 
the  surface  roughness  as  short  waves  become  steeper  on  the  swell  crest  than  those  on  the  swell  trough 
(Longuet-Higgins  and  Stewart,  1960).  The  purpose  of  the  proposed  study  is  to  quantify  the  spatial 
variability  of  the  ocean  surface  roughness  and  the  resultant  variations  of  momentum  and  energy 
exchanges  between  the  atmosphere  and  the  ocean  using  a  phase-resolving  numerical  wave  model. 

The  report  is  organized  as  follows.  First,  we  present  the  new  parameterization  of  the  momentum 
flux  transferred  from  the  wind  to  surface  gravity  waves  in  the  coastal  zone.  Next,  we  analyze  the 
effect  of  wave  nonlinearity  on  the  drag  coefficient  using  Stokes’  second-order  wave  theory.  After  that, 
we  show  howto  implement  the  parameterized  wind  stress  into  the  phase  resolving  Boussinesq  wave 
model  with  enhanced  dispersion  properties.  Then  tests  of  the  extended  Boussinesq  model  with  the 
wind  forcing  against  wave  measurements  in  a  shallow  lake  are  earned  out  followed  by  an  application 
of  the  model  to  the  study  of  sea-swell  interactions.  Finally,  we  summarize  the  findings  of  the  project. 


2.0  HETEROGENEOUS  AIR-SEA  MOMENTUM  FLUX 
2.1  Parameterization  of  Momentum  Flux 

The  momentum  flux  transferred  from  winds  to  surface  waves  in  the  time-domain  Boussinesq 
wave  model  may  be  parameterized  as  the  wind  stress,  T=paCdU \U\,  where  V  is  the  wind  velocity 

vector  at  a  reference  elevation,  and  Cd  is  the  corresponding  drag  coefficient.  Suggested  by  the  field 

evidence,  the  wind  drag  coefficient  in  the  coastal  zone  must  be  a  function  of  not  only  the  wind  speed, 
but  also  the  geometry  of  the  surface  gravity  waves  that  are  often  skewed  and  asymmetrical  owing  to 
the  seabed  boundary.  We  propose  a  new  formula  of  drag  coefficient  as  follows. 

Crf103  =  (a}  +a2  \  rjx  \)  +  bUl0  (1) 

where  a, ,  a2  and  b  are  empirical  coefficients,  and  7]x  is  the  instantaneous  surface  slope  computed  in 
the  Boussinesq  model.  Phase  average  of  the  drag  coefficient  allows  for  a  comparison  of  the  new 
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formula  with  conventional  drag  coefficient  formulas,  such  as  the  equation  proposed  by  Wu  (1980). 

Figure  1  illustrates  the  phase-averaged  drag  coefficient  as  a  function  of  wind  speed  and  wave 
steepness,  ka  (£= wave  number  and  a  =  wave  amplitude).  We  choose  a,  =0.2,  a2  =  18,  and  6=0.065. 
The  wave  steepness  ranges  from  0.05  to  0.2  with  an  increment  of  0.05.  It  is  seen  that  the  new  drag 
coefficient  increases  with  both  wind  speed  and  wave  steepness.  Following  Wu  (1980),  the  third  term 
in  Equation  (1)  is  set  to  be  0.4875  for  wind  speed  less  than  7.5m/s.  The  dashed  line  m  Figure  1 
depicts  the  drag  coefficient  given  by  Wu 

Crf103  =0.8  +  0.065I710  (2) 

which  is  only  a  function  of  wind  speed.  The  curves  of  the  new  drag  coefficient  are  parallel  to  Wu’s 
curve  and  the  offset  depends  on  the  wave  steepness. 

In  addition  to  Wu  (1 980)’  s  equation,  a  large  number  of  empirical  formulas  for  the  drag  coefficient 
have  been  proposed  in  the  literature  (see  e.g.  Geemaert,  1 990).  Each  formula  represents  the  best  fit  to 
a  specific  dataset.  Figure  2  shows  a  comparison  of  the  new  drag  coefficient  and  the  existing 
formulations.  It  is  seen  that  the  scattering  of  the  drag  coefficient  as  a  function  of  the  wind  speed  may 
be  explained  by  the  effect  of  the  sea  state,  or  the  wave  steepness.  Statistically,  Equation  (1)  agrees 
fairly  well  with  all  the  data,  which  is  not  achieved  by  any  single  existing  formula  of  drag  coefficient. 

One  of  the  crucial  differences  in  airflow  over  steep  surface  waves  in  comparison  to  the  flow  over 
a  flat  surface  is  the  presence  of  a  form  drag  due  to  the  separation  of  airflow.  Laboratoiy  measurements 
of  the  tangential  and  normal  stresses  on  the  air-water  interface  (Banner  and  Peirson,  1998)  suggest 
that  the  form  drag  accounts  for  the  major  proportion  of  the  wind  stress  once  the  waves  have 
developed  beyond  their  early  growth  stage.  The  most  interesting  scenario  of  wind  effects  on  coastal 
waves  is  the  strong  wind  condition.  We  hypothesize  that  the  separation  of  airflow  occurs  on  the  lee 
side  of  the  wave  crest  in  the  shallow  water  and,  consequently,  the  form  drag  is  a  dominant  contributor 
to  the  momentum  flux  between  the  atmosphere  and  the  nearshore  wave  field.  The  hypothesis  is  in  line 
with  the  sheltering  theory  proposed  by  Jeffreys  (1925)  as  a  mechanism  of  wave  generation. 

The  time-domain  Boussinesq  model  resolves  the  individual  wave  motion  and  predicts  reasonably 
well  the  geometry  of  the  surface  waves  distorted  by  the  seabed,  such  as  the  skewness  and  asymmetry 
of  the  waves  in  the  shallow  water  (e.g.  Chen  et  al„  2000  and  Shi  et  al.  2003).  The  wind  stress 
implemented  into  the  Boussinesq  model  should  vary  over  a  wavelength,  with  a  larger  drag  coefficient 
on  the  wave  crest  than  that  in  the  trough  to  take  the  form  drag  into  account.  In  fact,  only  the 
heterogeneous  distribution  of  the  wind  stress  over  a  wavelength  allows  for  changes  of  wave  height  by 
a  wind  in  a  phase-resolving  wave  model.  Assuming  that  the  form  drag  is  dominant  over  the  skm 
friction,  we  apply  the  wind  stress  on  the  wave  crest  only,  neglecting  the  effect  of  the  shear  stress  on 
the  wave  trough  in  the  Boussinesq  model. 

The  temporal  and  spatial  variation  of  drag  coefficient  is  one  of  the  important  features  of  the  wind 
stress  on  the  nearshore  waves  simulated  by  the  phase-resolving  Boussinesq  model.  With  the 
instantaneous  wave  celerity  and  wave  slope  obtained  from  the  computed  surface  elevation  at  each  aid 
and  every  time  step,  the  wind  stress  is  calculated  as  a  function  of  the  wind  velocity  relative  to  the 
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wave  celerity  as  well  as  the  wave  slope  in  the  following  form 

r  =  Cdpa  |  J/,o  -  C  |  (C/,o  -  C) 


(3) 


where  C  is  the  instantaneous  wave  celerity  estimated  in  the  Boussinesq  model.  Equation  0>)  is  similar 
to  the  wind  stress  formulation  followed  by  Schwab  et  al.  (1984)  in  a  wave  prediction  model  on  the 
basis  of  a  phase-averaged  momentum  balance  equation  rather  than  an  energy  balance  equation. 

2.2  Effects  of  Wave  Nonlinearity  on  the  Drag  Coefficient 

The  new  formulation  of  wind  drag  coefficient  introduced  in  the  preceding  section  depends  on  the 
wave  steepness  and  becomes  effective  only  on  the  wave  crest.  It  is  well  known  that  the  effects  of 
wave  nonlinearity  increase  the  wave  steepness  and  result  in  peaky  wave  profiles.  We  shall  analyze  the 
effect  of  wave  nonlinearity  on  the  drag  coefficient  given  by  Equation  (1)  using  Stokes’  second-order 

wave  theory. 

The  free  surface  profile  given  by  Stokes’  second-order  theory  reads 


^OCos(fa-^  +  — (2+C0Sh?)C°— C052(fa-mf)  (4) 

'  4  sinh  kh 

2  n  . 

where  a  is  the  first-order  (linear)  amplitude  of  the  free-surface  disturbance,  is  the  wave 

number,  h  is  the  still  water  depth,  L  is  the  wave  length,  x  is  the  horizontal  distance,  t  is  time,  and  co  is 
the  angular  frequency.  At  t  =  0  and  in  the  deep  water  ( h»L/2 ),  Equation  (4)  may  be  simplified  as 


n 


acoskx  + 


cos  2kx 


(5) 


Upon  the  use  of  the  trigonometric  identity  and  rearrangement  of  Equation  (5),  we  obtain 


r]k  =  a2k 2  cos2  kx  +  ak  cos  kx  - 


a2k2 

2 


(6) 


Because  the  wind  stress  (or  drag  coefficient)  is  only  applied  on  the  wave  crest  above  the  still  water 
level,  we  need  to  find  out  the  length  of  the  wave  crest.  Solving  Equation  (6)  with  rjk-0  yields 


x0k  =  arccos  / (ak) 


(7) 


in  which 


f(ak) =  - 


2  ak 


(— : )2+- 
y2akJ  2 


(8) 


Integrating  the  drag  coefficient  over  a  wave  length  leads  to  the  phase-averaged  diag  coefficient 
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c„ = i  f'c,& = + a.t,_  +  ^fa„c/,0  (?) 

d  L  *-^0  7C  n  It 

where  7;max  is  the  surface  elevation  at  the  wave  crest.  For  linear  waves,  we  have  x0=-  and  r]msx  =  a. 

The  ratio  of  the  drag  coefficient  on  a  Stokes’  second-order  wave  profile  to  that  on  a  linear  wave 
profile  is  given  by 

Cd  |  |  (^i +bf/10)(2/d'0 -7r)  +  flf2(^a)  (10) 

(Crf)to7~  +  (aI+bt/10)^  +  2a2to 

Figure  3  illustrates  the  variation  of  such  a  ratio  as  a  function  of  wave  steepness  and  wind  speed.  It  is 
seen  that  the  wave  nonlinearity  can  increase  the  phase-averaged  drag  coefficient  by  10-15  percent. 


AN  ADVANCED  BOUSS1NESQ  WAVE  MODEL  WITH  WIND  FORCING 

The  formulation  of  the  wind  stress  is  implemented  into  the  Boussinesq  wave  modeldeveloped  by 
Wei  et  al  (1995),  Kennedy  et  al.  (2000),  Chen  et  al  (2000),  and  Chen  et  al.  (2003).  In  order  to  take 
into  account  the  short  waves  generated  by  local  winds,  we  also  extend  thedispersion  accuracy  of  the 
one-dimensional  version  of  the  model  to  kh  =  6  (h  =  still  water  depth)  by  introducing  the  additional 
terms  from  Madsen  and  Schaffer’s  (1998)  equations  into  the  model.  A  new  set  of  coefficients  for 
better  shoaling  properties  of  the  equations  (Madsen,  2002,  personal  communication,  and  Kennedy  e 
al,  2002)  is  also  implemented  into  the  model.  In  one  horizontal  dimension,  the  extended  Boussmesq 

equations  read 


t],+Mx=0 


(11) 


where 


M  =  (h  +  rj ) 


ua  + 


fi  A(h2  -  + r  )k«, + {za  +  \(h -  w™*  L 
2  6  J  V  2  J 


(12) 


+ 


(&  -  7k  L +0*2- t*yn, -PifrnX 


and 


+  gVx  +A0  +A1  +A2  +A3  ~Rb  ~Rs+Rf  ~R «'  -0 


(13) 


where 


An  = 


^  +  {a2 -ax)h2  L*  +  (za  - a2h\hua )„  +  {a2  - a, )hhlxxx  - oc2h{h?lx)xx  (14) 
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A,  =  r.K„(K)„  +  f  W=,  -’?(*"«),  +  f  ((K)J 


A,  = 


^ (hua ),  -  T]Ua (hua L-^X 


^3  " 


1  2  (  \2  1  2 
2*7  (O  w«w« 


(15) 

(16) 

(17) 


In  the  equations,  r|  is  the  free  surface  elevation  relative  to  the  still  water  level;  ua  is  the  velocity  at  the 
reference  elevation  (=  -  0.54 1 22 h)  in  the  water  column;  and  the  subscripts  t  and  x  denote  time  and 
spatial  differentiations,  respectively.  The  dispersion  enhancement  coefficients  with  improved  shoaling 
properties  are(A,  A)  =  (°-03917’ °-315236)  and («,_«,)  =  (0.01052,  0.124537) .  The  additional 

terms,  Rb ,  Rs ,  Rf ,  and  Rw  represent  the  effects  of  wave  breaking,  subgrid  lateral  turbulent  mixing, 

seabed  shear  stress,  and  wind  stress,  respectively.  Detailed  descriptions  of  the  first  three  terms  can  be 
found  in  Chen  et  al.  (1999).  Upon  the  use  of  Equation  (3),  we  obtain 

i?i(,=-7^iCrf|t/10-C|([/10-C)  (18) 

p{h  +  J7) 

where  p  is  the  water  density. 

The  internal  wave-maker  for  the  generation  of  random  waves  at  the  offshore  boundary  is  extended 
to  accommodate  the  enhancement  of  dispersion  properties  in  the  Boussinesq  equations.  In  comparison 
to  Wei  et  al.’s  (1995)  model,  the  additional  terms  in  Equation  (12)  requires  the  solution  of  a 
tridiagonal  system  for  the  continuity  equation. 

Other  extensions  of  the  model  include  the  modification  of  the  breaking  criteria  in  the  model  to 
account  for  the  energy  dissipation  owing  to  whitecaps.  The  original  breaking  criteria  in  Kennedy  et  al. 
(2000)  were  designed  for  depth-limited  wave  breaking.  The  basic  idea  is  to  replace  the  still  water 
depth  h  by  CVg  where  g  is  the  gravitational  acceleration.  This  leads  to  the  modified  criteria  for  the 
onset  ( ?/,(/) )  and  cessation  (  rfp )  of  wave  breaking  as  follows 

=7]C  (19) 

V \F)=Yic 

in  which  y,  =0.35-0.65  and  y2  =0.1 5.  We  choose  the  lower  limit  of  y,  in  the  presence  of  a  wind  that 
is  in  the  direction  of  the  wave  train.  Other  terms  with  h  in  the  breaking  scheme  of  Kennedy  et  al. 
(2000)  are  also  modified  in  a  similar  fashion. 
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4.0  MODELING  WAVE  GROWTH  ON  A  SHALLOW  LAKE 

The  test  cases  in  Chen  et  al.  (2002)  were  focused  on  the  surface  roughness  length  and  drag 
coefficient  associated  with  shoaling  waves  over  a  distance  of  about  500m.  Can  the  wind  stress 
implemented  into  the  Boussinesq  model  simulate  the  growth  of  wind  waves  over  a  longei  fete  .  e 
shall  address  such  a  question  using  the  data  set  of  wave  growth  in  Lake  George,  Australia  collected  by 
Young  and  Verhagen  (1996). 

Lake  George  is  a  shallow  lake  with  a  typical  water  depth  of  2m.  It  is  about  20km  long  and  10km 
wide  A  series  of  eight  observation  stations  were  deployed  along  the  north-south  fetch  to  measure  the 
wind  waves.  The  wave  data  collected  by  Young  and  Verhagen  (1996)  under  a  nearly  ideal  condition 
have  served  as  a  test  bed  for  a  number  of  phase-averaged  wave  models,  such  as  SWAN  (Simulating 
WAves  Nearshore,  Booij  et  al.,  1999).  We  choose  two  data  sets  corresponding  to  a  medium  wind 
speed  of  U]0  =  10.8m/s  and  a  large  wind  speed  of  Uw  =  15.2m/s  to  test  the  wind  stress  formulation 

incorporated  into  the  Boussinesq  model.  Owing  to  the  dispersion  limit  of  the  Boussinesq  model  <kh  < 
6)  only  the  observations  at  the  last  four  stations  (5-8)  are  utilized.  The  fetch  between  Stations  P  and  8 
is  about  4.25  km.  The  up-wave  boundary  condition  is  taken  from  the  observation  at  Station  5  where 
the  zero-moment  wave  height  is  Hmo  =  0.367m  and  the  peak  wave  period  isTp  =  2.34s  in  the  case  of  a 

moderate  wind  (t/10  =  10.8m/s),  and  Hm  =  0.473m  and  T  =  2.4s  in  the  case  of  a  strong  wind  (N10  - 

1 5  2m/s)  We  use  TMA  shallow  water  wave  spectra  (Bouws  et  al.,  1985)  with  the  shape  parameter  y- 
3 .3  at  the  up-wave  boundary  as  the  input  to  the  Boussinesq  model.  A  bottom  friction  coefficient  of/- 
0.0005  is  used  in  the  quadratic  law  of  the  bottom  shear  stress  (Chen  et  al.,  1999). 

Figure  4  shows  comparisons  of  the  computed  and  measured  significant  wave  heights  as  well  as 
the  computed  and  inferred  drag  coefficients  under  the  moderate  wind  condition.  It  is  seen  that  the 
wave  heights  predicted  by  the  Boussinesq  model  are  in  fairly  good  agreement  with  the  field 
measurements.  The  computed  peak  wave  periods,  however,  do  not  agree  with  the  observations.  The 
wind  stress  incorporated  into  the  Boussinesq  model  simply  does  not  lead  to  the  increase  of  wave 
period  along  the  fetch.  The  downshift  of  the  peak  frequency  is  absent  in  the  modeled  wave  energy 
spectrum.  This  is  attributed  to  the  inability  of  the  on^dimensional  Boussmesq  model  to  take 
quadruplet  interactions  into  account.  Interestingly,  the  phase-averaged  drag  coefficient  calculated 
from  the  Boussinesq  model  is  in  good  agreement  with  the  coefficient  estimated  using  themo  e  e 
wave  characteristics  and  Anctil  and  Donelan’s  (1 996)  formula.  The  drag  coefficients  given  by  Taylor 
and  Yelland’s  (2001)  formula  (triangles),  however,  are  smaller  than  either  Boussmesq  or  Anctil  and 

Donelan’s  result  (solid  line  and  squares). 

Similar  comparisons  are  made  for  the  case  of  large  wind  speed  as  shown  in  Figure  5.  Again,  the 
Boussinesq  model  predicts  the  growth  of  wave  height  fairly  well  in  comparison  to  the  field 
measurements  but  considerably  under-predicts  the  increase  of  wave  penod.  The  agreement  among  the 
drag  coefficients  is  also  very  similar  to  the  case  with  a  moderate  wind.  The  phase-averaged  drag 
coefficient  obtained  from  the  Boussinesq  model  agrees  better  with  Anctil  and  Donelan’s  result  than 
does  Taylor  and  Yelland’s  formula,  as  shown  in  Figure  5b.  Notice  that  all  the  drag  coefficients  are 
computed  using  the  wave  field  given  by  the  one-dimensional  Boussinesq  model,  which  is  not  able  to 
predict  the  downshift  of  the  peak  frequency.  Nevertheless,  the  prediction  of  wave  height  growth  under 
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both  moderate  and  strong  wind  conditions  by  the  Boussines,  mode)  confirms  that  the  new 
formulation  of  wind  stress  is  a  good  representation  of  the  momentum  flux  transferred  from  a  wmd  to 
surface  waves  in  a  phase-resolving  model. 


5.0  THE  INFLUENCE  OF  SWELL  ON  THE  GROWTH  OF  WIND  WAVES 

Donelan  (1987)  among  others  has  found  that  the  addition  of  mechanically-generated  swell  to  the 
wind  waves  in  a  wave  flume  results  in  a  pronounced  reduction  in  the  energy  of  the  wmd  sea  while  the 
swell  gains  significant  energy  from  the  following  wind  in  comparison  to  pure  wind  waves^  Although  a 
number  of  hypotheses  have  been  proposed  in  the  literature,  the  attenuation  mechanism  for  the  wind 
sea  owing  to  the  presence  of  swell  is  still  poorly  understood.  As  demonstrated  m  the  preceding 
section,  the  extended  Boussinesq  model  is  able  to  simulate  the  growth  of  wave  energy  for  a  given 
wind  speed  and  a  non-zero  wave  spectrum  at  the  up-wave  boundary.  We  shall  conduct  a  senes  of 
numerical  experiments  to  examine  how  swell  influences  the  growth  rate  of  wind  waves  on  the  basis  of 
the  extended  Boussinesq  wave  model  with  the  wind  forcing.  The  Lake  George  case  with  a  moderate 
wind  is  used  as  a  reference  for  comparisons. 

Fiaure  6  depicts  snapshots  of  the  computed  free  surface  elevations  of  the  cases  with  and  without 
swell°The  right  panel  corresponds  to  the  first  case  with  a  moderate  wind  speed  m  the  preceding 
section.  A  monochromatic,  long  wave  train  with  a  period  of  three  times  the  peak  wave  penod  of  the 
short  waves  is  imposed  at  the  up-wave  boundaiy  in  addition  to  the  wmd  waves,  as  shown  m  the  right 
panel.  The  wave  height  of  the  swell  is  0.184  m.  Each  snapshot  (e.g.  A  to  F)  m  each  panel  is  15 
minutes  apart,  beginning  15  minutes  after  the  model  starts. 

Fiaure  7  shows  the  energy  density  spectra  of  the  modeled  sea  surface  with  and  without  swell.  The 
dashed  lines  represent  the  spectra  at  the  up- wave  boundary  while  the  solid  lines  are  the  spectra  4 
down-wind.  In  contrast  to  the  spectra  of  the  pure  wind  waves  on  the  right,  the  left  panel  shows  that 
the  short  waves  barely  grow  while  the  swell  gains  energy  from  the  wmd.  The  Boussinesq  model 
qualitatively  reproduces  what  Donelan  (1987)  observed  in  the  laboratory  about  the  suppression  of  the 
wave  growth  rate  in  the  presence  swell  under  field  conditions.  No  field  measurements  of  sea-swell 
interaction,  however,  are  available  for  a  direct  test  of  the  influence  of  swell  predicted  by  the 
Boussinesq  wave  model  with  the  wind  forcing. 

The  slope  of  the  sea  surface  is  a  measure  of  the  surface  roughness.  Figure  8  illustrates  spectra  of 
the  modeled  sea  surface  slope  with  and  without  swell.  The  dashed  lines  represent  the  slope  spectra  at 
the  up-wave  boundary  while  the  solid  lines  are  the  spectra  4  km  down- wind.  A  comparison  of  the  left 
and  right  panels  shows  that  the  swell  actually  reduces  the  growth  rate  of  the  surface  slope,  or  the  sea 
surface  roughness  of  the  short  waves.  This  is  consistent  with  the  suppression  of  the  shortwave  energy 
growth  rate  by  the  following  swell. 

Figure  9  depicts  the  spatial  variations  of  the  wave  height  and  phase-averaged  drag  coefficient  with 
and  without  swell.  The  solid  lines  denote  the  model  result  of  pure  wind  waves  while  the  dashed  lines 
are  the  result  of  combined  sea  and  swell.  It  is  seen  that  the  presence  of  swell  yields  a  slightly  smaller 
growth  rate  of  the  total  wave  energy  than  that  of  the  pure  sea.  A  similar  trend  is  found  in  the  phase- 
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averaged  drag  coefficient,  as  shown  on  the  right.  To  examine  the  effect  of  the  swell  steepness  on  the 
growth  rate  of  the  wind  waves,  we  reduce  the  swell  amplitude  to  0.041  m,  and  repeat  the  simulation 
above  without  changing  other  conditions.  The  Boussinesq  model  predicts  that  the  smaller  the  swell 
steepness,  the  less  degree  the  suppression  of  the  shore-wave  growth  rate,  as  seen  by  a  direct 
comparison  of  Figure  7  and  Figure  10. 

The  numerical  results  suggest  that  on  the  swell  crests,  the  larger  surface  roughness  and  the 
sheltering  effect  lead  to  a  larger  momentum  transfer  than  that  on  the  swell  trough,  but  the 
momentum/energy  flux  averaged  over  the  swell  wavelength  remains  relatively  unchanged.  The 
majority  of  the  short  waves  riding  on  the  swell  crest  are  subjected  to  the  swell-induced  form  drag 
owin'*  to  the  sheltering  effect  and  a  homogeneous  wind  stress  over  an  entire  short  wavelengt 
prevents  the  short  wave  from  growing.  Therefore,  the  swell  takes  advantage  of  the  stronger 
momentum  flux  on  its  crest  and  harvests  most  of  the  energy  transferred  from  the  wind  to  the  water.  In 
layman’s  words,  the  finite  amplitude  swell  serves  as  a  “master”  wave  and  the  wind  waves  riding  on 
the  swell  act  as  “slaves”  for  momentum  transfer.  Consequently,  the  swell  amplitude  increases  while 

the  wind  sea  barely  grows. 


6.0  SUMMARY  AND  CONCLUSIONS 

A  phase-resolving  wind  stress  formula  has  been  developed,  analyzed  and  applied.  The  report 
documents  (1)  the  parameterization  of  the  momentum  flux  transferred  from  the  wind  to  surface 
gravity  waves  in  a  phase-resolving  Boussinesq  wave  model,  (2)  the  analysis  of  the  effect  of  wave 
nonlinearity  on  the  drag  coefficient  based  on  Stokes’  second-order  wave  theory,  (3)  tests  of  the 
extended  Boussinesq  model  with  the  wind  forcing  against  wave  measurements  in  a  shallow  lake,  and 

(4)  numerical  experiments  on  the  influence  of  swell  on  the  growth  rate  of  the  wind  waves.  Fairly  good 

agreement  between  the  modeled  wave  heights  and  field  data  has  been  found.  A  new,  simple 
mechanism  has  been  suggested  to  explain  the  suppression  of  the  wind  wave  growth  rate  in  the 
presence  of  swell.  Results  of  the  study  have  been  published  in  the  J  oumal  of  Waterway ,  Port,  coastal 
and  Ocean  Engineering  (Chen  et  al.  2004a)  and  the  International  Conference  on  Coastal  Engineering 
(Chen  et  al.  2002  and  Chen  et  al.  2004b). 

The  methodology  for  the  parameterization  of  the  phase-resolving  air-sea  momentum  flux  as  well 
as  the  extended  Boussinesq  wave  model  incorporating  the  wind  effects  appear  to  be  a  promising  tool 
for  the  study  of  sea-swell-wind  interactions  and  wind  effects  on  nearshore  wave  propagation 
However,  the  limit  of  the  dispersion  accuracy  to  kh  <  6  prevents  the  extended  Boussinesq  wave  mo  e 
from  being  compared  directly  with  existing  laboratory  datasets  on  the  influence  of  swell  on  the 
growth  of  wind  waves.  Thus,  it  is  recommended  that  state-of-the-art  Boussinesq  models  with  the 
dispersion  accuracy  up  kh  =  40  (e.g.  Madsen  et  al.,  2002)  be  employed  for  the  study  of  sea-swell-wind 
interactions  in  the  future. 
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9.0  APPENDIX 


Fig.  1  Drag  coefficient  as  a  function  of  wind  speed  and  wave  steepness  (solid  lines).  The  dashed 
line  is  Wu’s  (1 980)  drag  coefficient. 


14 


2  Comparison  of  Equation  (1)  with  the  drag  coefficient  formulas  complied  by  Geemaert  (1 990). 


L5 


(a)  Comparison  of  Wave  Height 


Distance  (km) 


FK  4  Comparisons  of  the  computed  and  measured  significant  wave  heights  (top),  and  the  computed 
and  inferred  drag  coefficients  (bottom)  in  the  case  of  Ul0=  10.8m/s.  Solid  lines:  Boussinesq  model 
results,  circles:  Young  and  Verhagen’s  (1 996)  observations,  squares:  Anctil  and  Donelan  (1 996),  and 
triangles:  Taylor  and  Yelland  (2001). 
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(a)  Comparison  of  Wave  Height 


Distance  (km) 


FK  5  Comparisons  of  the  computed  and  measured  significant  wave  heights  (a)  as  well  as  the 
computed  and  inferred  drag  coefficients  (b)  in  the  case  of  Uw=  15.2m/s.  Solid  lines:  Boussmesq 
model  results,  circles:  Young  and  Verhagen’s  (1996)  observations,  squares:  Anctil  and  Donelan 
(1996),  and  triangles:  Taylor  and  Yelland  (2001). 
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Wind  Waves  without  Swell 


Wind  Waves  with  Swell 


Fig.  6  Snapshots  of  modeled  sea  surface  without  (left)  and  with  (right)  swell 
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Fig.  7  Suppression  of  wave  growth  by  the  ambient  swell  (left).  Dashed  lines:  up-wave  boundary,  solid 
lines:  4  Ion  down- wave. 
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(b)  Wind  Waves  without  Swell 


Frequency  (Hz) 


Fig.  8  Spectra  of  sea  surface  slope  with  (left)  and  without  (right)  swell.  Dashed  lines,  up-wave 
boundary,  solid  lines:  4  km  down-wave. 


21 


Hs(m> 


Comparison  of  Wave  Height 


Fig.  9  Effects  of  swell  on  wave  height  and  drag  coefficient.  Dashed  lines:  sea  and  swell,  solid  lines: 
short  waves. 


Fig.  10  Suppression  of  wave  growth  by  the  swell  (left)  with  a  small  amplitude.  Dashed  lines,  up-wave 
boundary,  solid  lines:  4  km  down-wave. 


23 


